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We consider the construction of nonsingular Pre-Big-Bang and Ekpyrotic type cosmological 
models realized by the addition to the action of specific higher-order terms stemming from quantum 
corrections. We study models involving general relativity coupled to a single scalar field with a 
potential motivated by the Ekpyrotic scenario. We find that the inclusion of the string loop and 
quantum correction terms in the string frame makes it possible to obtain solutions of the variational 
equations which are nonsingular and bouncing in the Einstein frame, even when a negative expo- 
nential potential is present, as is the case in the Ekpyrotic scenario. This allows us to discuss the 
evolution of cosmological perturbations without the need to invoke matching conditions between two 
Einstein Universes, one representing the contracting branch, the second the expanding branch. We 
analyze the spectra of perturbations produced during the bouncing phase and find that the spectrum 
of curvature fluctuations in the model proposed originally to implement the Ekpyrotic scenario has 
a large blue tilt (n-R, — 3). Except for instabilities introduced on small scales, the result agrees 
with what is obtained by imposing continuity of the induced metric and of the extrinsic curvature 
across a constant scalar field (up to k^ corrections equal to the constant energy density) matching 
surface between the contracting and the expanding Einstein Universes. We also discuss nonsingular 
cosmological solutions obtained when a Gauss-Bonnet term with a coefficient suitably dependent on 
the scalar matter field is added to the action in the Einstein frame with a potential for the scalar 
field present. In this scenario, nonsingular solutions are found which start in an asymptotically fiat 
state, undergo a period of super-exponential inflation and end with a graceful exit. The spectrum 
of fluctuations is also calculated in this case. 



I. INTRODUCTION 

There has recently been a lot of interest in cosmological scenarios in which it is assumed that instead of emerging 
from an initial Big Bang singularity, our Universe has resulted from an Einstein frame bounce which connects a 
previous contracting phase with the present phase of cosmological expansion. A lot of this interest has been fueled 
by string cosmology, the attempt to merge string theory and cosmology. Pre-Big-Bang (PBB) cosmology [1,2] (see 
[3,4] for a comprehensive review) and the Ekpyrotic scenario [5] are two well-known models in which our present 
phase of cosmological expansion is postulated to have emerged from a previous phase of cosmological contraction ^. 
In both examples, however, the cosmological description in terms of an effective action breaks down at the bounce. In 
the case of PBB cosmology this bounce corresponds to a region of high curvature where higher derivative and string 
corrections to the effective action will be important, in the case of the Ekpyrotic scenario the bounce occurs when 
two four space-time dimensional branes collide in a five dimensional bulk. 



^In PBB cosmology this statement is true from the point of view of the Einstein frame metric, in the Ekpyrotic scenario it is 
true from the point of view of the four space-time dimensional effective action which is used to describe the cosmology. 
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Models with a cosmological bounce potentially provide an alternative to cosmological inflation in addressing the 
homogeneity problem of standard cosmology and in yielding a causal mechanism of structure formation, the latter 
since at times long before the bounce fixed comoving scales of cosmological interest today will have been inside 
the Hubble radius.^ However, since in both PBB and Ekpyrotic scenarios the Hubble parameter increases during the 
collapsing phase, symmetry arguments such as those used originally [12] to predict the scale- invariance of cosmological 
fluctuations in inflationary cosmology would lead one to expect a blue spectrum of curvature perturbations in these 
models, at least in effective field theory models in which there is only one "matter" field. As outlined in Appendix A, 
in PBB cosmology one expects a spectrum with spectral index n = 4, whereas in the Ekpyrotic scenario one expects 
n = 3. In the case of PBB cosmology, this heuristic prediction was confirmed [13] by a general relativistic analysis 
(which is, however, subject to the caveats indicated below). In the case of the Ekpyrotic scenario, there is a large 
disagreement in the results. Whereas the works of [14-17] yield results in agreement with the heuristic prediction 
(namely n = 3), others [18,19] obtain a scale-invariant spectrum of adiabatic fluctuations (thus also casting doubt on 
past results in the literature on the spectrum of fluctuations in the PBB scenario). 

The singularity of the effective action at the time of the bounce makes it impossible to follow the evolution of the 
background cosmology and of the resulting cosmological perturbations rigorously [20] . In much of the previous work, 
both in the context of PBB cosmology [21] and of the Ekpyrotic scenario [15-17,19], the fluctuations were computed 
by matching two Einstein Universes (the first representing the contracting phase, the second the expanding phase) 
along a space-like surface (representing the bounce region) and applying continuity of the induced metric and of the 
extrinsic curvature across the surface [22,21]. As emphasized in [19], the result will depend on how the matching 
surface is chosen ^. 

In the context of PBB cosmology, it was realized [26 28] (see also [29-31]) that higher derivative corrections 
(defined in the string frame) to the action induced by inverse string tension and coupling constant corrections can 
yield a nonsingular background cosmology This then allows the study of the evolution of cosmological perturbations 
without having to use ad hoc matching prescriptions. The effects of the higher derivative terms in the action on the 
evolution of fluctuations in the PBB cosmology was investigated in [36]. It was found that for low frequency modes, 
the spectrum of fluctuations is unaffected by the higher derivative terms, and the result obtained is the same as what 
follows from the analysis using matching conditions between two Einstein Universes [13,21] joined along a constant 
scalar field hypcrsurfacc. 

Since the Ekpyrotic scenario makes use of a negative exponential potential for the scalar matter field, which leads 
to an extra instability of the system, it is not clear that the higher derivative terms used in [26-28] can in this case 
achieve a nonsingular cosmology. The first main result of this paper is that, with suitably chosen coefficients, the 
above mentioned terms are indeed sufficient to produce a nonsingular cosmology. 

In this paper, we add the same higher derivative terms used in [27] to the action (which includes a positive or negative 
potential for the scalar matter field) and construct nonsingular bouncing cosmologies. At the level of the effective 
action, our Lagrangian can be viewed as giving both nonsingular solutions of modified PBB type (the modification 
consisting of the addition of an exponential potential for the dilaton), and also nonsingular Ekpyrotic solutions. The 
justification for adding these higher derivative terms is different in the cases of modified PBB cosmology and in the 
Ekpyrotic scenario. In the case of PBB cosmology, both the string coupling constant and the curvature become large 
as the dilaton increases, thus justifying the inclusion of both higher derivative terms of the gravitational action and of 



^See refs. [611] for critical arguments on the Ekpyrotic scenario. 

^As emphasized already in [15] and [23], there is a consistency check for proposed matching surfaces: when applied to the 
reheating surface in inflationary cosmology, the correct result should emerge. This does not happen with the prescription 
advocated in [19] , nor does it with the matching prescription of [18] which is not based on a geometric analysis (see also [24,25] 
for a criticism of the latter matching proscription). 

■^Construction of nonsingular cosmologies in pure Einstein gravity by means of specific higher derivative terms is also possible 
(see e.g. [32,33] and its application to PBB cosmology in [34,35]). 
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quantum corrections. In the case of Ekpyrotic cosmology (wc have the initial scenario of [5] in mind in which a bulk 
brane impacts our physical space-time orbifold fixed plane at the time of the bounce and in which the dilaton and 
hence the string coupling constant are fixed), the density and hence curvature at the bounce are large, thus justifying 
including higher derivative terms. In addition, the brane collision is a quantum mechanical process, thus justifying 
including loop corrections in the action. Note that our method yields a way of constructing a nonsingular bouncing 
Universe which works even in a spatially flat Universe and is thus different from the constructions of [37] in which a 
positive spatial curvature is used to generate a bouncing cosmology. Since tracing back the spatial curvature into the 
very early Universe given the present date leads - under the assumption that there was no period of inflation after the 
bounce - to a highly suppressed curvature at early times, our approach in obtaining a bouncing cosmology appears 
more realistic. 

We follow the fluctuations through the bounce, and study the spectrum of the resulting cosmological perturbations 
at late times. In this analysis, no matching conditions at the bounce are necessary. Note, however, that in principle 
the final spectrum could depend on the frame in which the higher order correction terms are introduced, and on the 
specific form of the correction terms. In our nonsingular scenario discussed in Sec. IV, the correction terms are defined 
in the string frame and we find that the final spectrum of cosmological fluctuations on long wavelength scales has 
a shape which agrees with what is obtained when applying the matching conditions of [22,21] on a constant scalar 
field surface ^ (the most physical choice of the matching surface both in PBB and Ekpyrotic models). In particular, 
for the Ekpyrotic model of [5] rendered nonsingular by our construction, we obtain a blue spectrum of the curvature 
perturbation with index n-ji = 3. 

II. SINGLE SCALAR FIELD WITH AN EXPONENTIAL POTENTIAL 

The Lagrangian considered in this paper can be used to describe both a modified PBB model in which the dilaton 
has an exponential potential as well as the Ekpyrotic scenario. Our Lagrangian describes gravity plus a single scalar 
matter field (f>. In the case of PBB cosmology, the physical frame is the string frame, and (f> is the dilaton field. In the 
case of the original version of the Ekpyrotic scenario [5] , the physical frame is the Einstein frame since the dilaton is 
fixed, and the field ^ is related to the separation of a bulk brane from our four-dimensional space-time orbifold fixed 
plane. In the case of the second version of the Ekpyrotic scenario [38] and in the cyclic variant thereof [39,40], ^ is 
the modulus field denoting the size of the orbifold (the separation of the two orbifold fixed planes). 

We begin with the Lagrangian of the four-dimensional effective theory in the string frame, which is: 



lR+l{yct>f-Vs{<P) 



(2.1) 



where R is the Ricci scalar, Vs{(t>) is the scalar field potential in the string frame. In this form, the action looks 
reminiscent of the action for PBB cosmology. Note that Vs{4') = in the simplest version of the PBB scenario. 
We set the units such that SttG = 1 with G being a four-dimensional gravitational constant. Making a conformal 
transformation 

g^u = e-'I'g^. , (2.2) 
the action in the Einstein frame can be written as 

'1^1 



Se = I d x^-g 



2i?--(V</))^-14;(</.) 



(2.3) 



^Note that up to terms of order fe^, the constant scalar field surface and the constant energy density surface are identical, as 
discussed in [15]. 
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where 



Introducing a rescaled field = ±(j)/V2, the action (2.3) reads 



(2.4) 



(2.5) 



In this form, the action is seen to describe both the PBB model in the Einstein frame, as well as the Ekpyrotic scenario 
[19]. 

The Ekpyrotic scenario is characterized by an exponential potential [5] 



Ve = -Vo exp 



P 



(2.6) 



with < p <^ 1. The field (p denotes the separation of two parallel branes. According to the Ekpyrotic scenario, 
the branes are initially widely separated but are approaching each other, which means that ip begins near +oo and 
is decreasing toward (p = 0. In the PBB scenario, in contrast, the dilaton starts out from a weakly coupled regime 
with (f) increasing from — oo. Thus, if we want the potential (2.6) to describe a modified PBB scenario with a dilaton 
potential which is important when </> — > but negligible for — > — oo, we have to use the relation (p = — 0/\/2 between 
the field p in the Ekpyrotic case and the dilaton (j> in the PBB case. 

Adopting the Friedmann-Robertson-Lemaitre- Walker (FRLW) metric ds"^ = —dt% + a%dx% in the Einstein frame, 
the background equations are given by 



3HI = ^cp' + Ve{^), <p + 3He^ + V^{^)=0, 



(2.7) 



where a prime denotes a derivative with respect to a cosmic time, ts- For the exponential potential (2.4) we have the 
following exact solution 



aE oc [-tsr, He = — , Ke = 75 

lE ti 



f2 



2p 



Ie 



(2.8) 



The solution for < describes the contracting universe in the Einstein frame prior to the collision of branes. 
In the string frame the action is given by (2.1) with potential 



Vs = -Vbexp 



The FRLW metric in the string frame is described by ds^ 
in the Einstein frame as 



1 

VP 



(2.9) 



-dt'g + a.gdx'g), which is connected to the quantities 



dts 



e-'^'^dtE, as = e-'Pl^aE , 



where we used the relation (p = —\f2ip>. Integrating the first relation gives 

-{\-^)ts^{-tEf-^ . 

Therefore the evolution of 05 and ^ in the string frame is 



as 



ln[-(l-Vp)is] 



(2.10) 



(2.11) 



(2.12) 



This illustrates the super- infiationary solution with growing dilaton from (j) = —00. Note that singularities are 
inevitable in both frames as t — > 0. We wish to analyze whether this singularity can be avoided by including higher- 
order corrections. 
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III. GENERAL ACTIONS AND EVOLUTION EQUATIONS 



In this section we present the background and perturbed equations in the case of a generahzed action containing 
higher derivative terms. We write this action in the form [36,41] 



-J 



i/(i?,(/.)-^a;(<^)(V<^)2-y(.^)+£, 



(3.1) 



where f{R, cj)) is a function of the Ricci scalar R and a scalar field (j). oj{(j}) and V{(j)) are general functions of (j). The 
Lagrangian Cc represents the higher-order corrections to the tree-level action. Both higher derivative gravitational 
terms and terms involving Lp appear. The action (3.1) applies not only to low-energy effective string theories, but also 
to effective action approaches to Einstein quantum gravity and to scalar tensor theories, among others. 

As mentioned in the Introduction, our motivation for considering the addition of higher derivative terms in the 
effective action is to construct a nonsingular bouncing model and thus to overcome the singularity problem ( "graceful 
exit problem") which plagues both the PBB and the Ekpyrotic scenario. The higher-order contribution Cc can be 
written as the sum of the a' classical correction and the quantum loop correction Cq [27,42]. Both involve the 
same gravitational and scalar field terms, but are multiplied by different powers of e~'^. 

The leading a' (string) correction to the gravitational action we adopt is given by [26,27] 

= - \a'\m [cRl^ + d{ycj>f] , (3.2) 

where £^{<f)) is a general function of (j) and iJ^g = R? — AR'^" R^^ + R^"°'^ Rixvai3 is the Gauss-Bonnet term. The inverse 
string tension, a', is set to unity. The Gauss-Bonnet term has the property of keeping the order of the gravitational 
equations of motion unchanged. It has been known since the early days of string theory that this term arises as the 
lowest string correction to the gravitational field equations in a string theory background. At tree level (lowest order 
in h), we have ^((/)) = —e~'^. When applied to PBB cosmology, it turns out that as t ^ — oo, the invariants Rq^ 
and (V(?i)^ decay faster than the coefficient function ^{(p) blows up. Hence, the correction terms in the Lagrangian 
are unimportant for large negative values of (p, but become important as the system approaches the strongly coupled 
region (^ ~ 0). 

Following [27] , we take the higher n-loop correction terms Cq in addition to the tree-level term . For the moment, 
however, we will keep ^{(f)) general. We will give specific forms for ^{<p) and Cq later. 



A. Background equations 

Variation of the action (3.1) with respect to the scale factor, the lapse function (then set to 1 after the variation) 
and the scalar matter field leads to the following background equations [36] 

H^ = -^{oj^^+RF-f + 2V-6HF + p,), (3.3) 

^ = i f-^'^' + HF-F-Ipc- Ipc) , (3.4) 



2F V ^ 2^ 2' 



4> + 3H^+^(uJ^<j>^-f4, + 2V4,-A4,'^=0, (3.5) 

where F = df/dR and H = a/a. pc, Pc, and correspond to the higher-order curvature and derivative corrections 
with stress-energy tensor Tjf = {—PcPcPcPc)- A^ comes from the variation of Cc with respect to <p. For the tree-level 
a' correction (3.2), one has 

Pc = 2a' X (l2ciH^ - \di4>''^ , (3.6) 
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Pc = -2a' X { 4c 
Arf, = -a'X 



iH' + 2^H{H + H') + } , (3.7) 



2 

2ici^H'\H + H'^) - d(f\3^(p + 12^(j) + U^^H) . (3.8) 



Note that taking into account quantum loop corrections provides additional source terms for pc, Pc, and A^. We will 
discuss this issue in Sec. IV. 

B. Perturbation equations 

A perturbed space- time metric has the following form for scalar perturbations in an arbitrary gauge (see e.g. [43], 
where the function A is denoted by (p): 

ds"^ = -(1 + 2A)dt^ + 2a{t)B^idx'dt + a'^{t)[{l - 2ip)Sij + 2E^ij]dx'dx^ , (3.9) 

where a comma denotes the usual flat space coordinate derivative. We introduce the curvature perturbation, TZ, in 
the comoving gauge [44] 

H = ii) + ^6(j). (3.10) 
The perturbed Einstein equations for the action (3.1) are written in the form [41,36] 

J^(a3Q7t)--47e = 0, (3.11) 

where 

{H + lf ' ^ ' ' 

^ ^ ^ ^ 4AcC.j>^ - IQXcjHI + 8Ac(| - iH)P 
uj^'^ + 3I{F-4:X(4H'^)-6Xd^4)^ 

with I={F- 4XciH^)/{2F - SAc^iJ). 

Introducing a new quantity, ^ ^ = zTZ, with z = a\/Q, each Fourier component of ^ satisfies the second order 
differential equation 

*ft+(sfc'-^)*fc=0> (3.14) 

where a prime denotes the derivative with respect to conformal time, V = J a~^dt. In the large scale limit, jsfc^j <C 
\z"/z\, eq. (3.14) is integrated to give 



T^k = Ck + Dk 



(3.15) 



where Ck and Dk are integration constants. The curvature perturbation is conserved on super-Hubble scales as long 
as the second term in eq. (3.15) is not strongly dominating, as in the case of the single field, slow-roll inflationary 
scenarios. 



Note that if is the variable in terms of which - for unmodified Einstein gravity - the action for fluctuations has the canonical 
form of a free field action with time dependent mass (see e.g. [43] where the variable is denoted as v). 
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If the evolution of z before the bounce is given in the form 



z^{-nr, (3.16) 

the second term in eq. (3.15) yields ^ drj/z^ oc (— ??)^~^^. Therefore curvature perturbations can be amplified for 
7 > 1/2 on super-Hubble scales, while they are not for 7 < 1/2 [17] (Note that TZk oc ln(— 77) for 7 = 1/2). Whether 
this enhancement occurs or not depends on the time evolution of z, and therefore on the string cosmological model. 

We need to go to the next order solution of oq. (3.14) in order to obtain the spectrum of curvature perturbations. 
If s is a positive constant (as it will be in the asymptotic limits) , the solution for is expressed by the combination 
of the Hankel functions: 



where 



c^Hi^\x) + C2H'^\x)\ , (3.17) 



x=v^A:|r;| v=^\l-2^\. (3.18) 



The solution (3.17) corresponds to the Minkowski vacuum state in the small-scale limit {k 00). 
We can expand the Hankel functions in the following form [45]: 



" ^n!V 4/ siuTTi/ \2J T(iy + n+l) \2J T{-v + n- 



n=0 



(3.19) 



The curvature perturbation, TZk = "^k/z, has two solutions which are proportional to k''\rj\''~'^~^'^/'^ and 
k~'^\ri\~^~^'^'^/'^ , which follow from the first and second term in eq. (3.19), respectively. In the large-scale limit 

{k 0), the contribution of the second term dominates over the first one, thereby yielding the spectrum of the 
curvature perturbation as 

= 1^ iT^fcl' oc /c3-2- oc (3.20) 

in which case the spectral tilt is 

riK- 1 =3- |1 - 27I . (3.21) 

Note that we have k~''\r}\~''~''~^^/'^ = k~''\r]\^ for 7 < 1/2, in which case the constant mode Ck in eq. (3.15) corresponds 
to the solution which comes from the second term in eq. (3.19). The Ekpyrotic scenario with a negative potential (0 < 
p < 1/3) belongs to this case (7 < 1/2) as we will show later. When 7 > 1/2 one has k-''\r]\-''-^+'^/'^ = Aj-'^ItjI^-^t, 
which means that Tlk grows before the graceful exit (77 < 0). We will discuss a string- inspired model that belongs to 
this case in Sec. V. The PBB scenario corresponds to the marginal case with 7 = 1/2. 

Note that s is exactly unity in eq. (3.14) when the corrections Lc are not taken into account. In the presence of 
higher-order corrections [Cc ^ 0), s is generally a time-varying function, in which case the formula (3.21) can not be 
directly applied. Nevertheless it is still valid if s is a slowly varying positive function. 

In subsequent sections we shall apply the above general formulas to concrete string-inspired models. In Section 
IV we apply the string loop and quantum corrections to the low energy effective action in the string frame and 
find nonsingular bouncing cosmological solutions. In Section V, we consider a situation with fixed dilaton and add 
a Gauss-Bonnet term to the Einstein frame action. We find nonsingular cosmological solutions which begin in an 
asymptotically flat state, undergo a period of super-exponential inflation which terminates with a graceful exit. 



'^Note that as long as the additional terms Cc in the action are negligible, then 7 = p/(l — p), and thus < 7 <C 1 for the 
collapsing phase of Ekpyrotic cosmology, and 7 = 1/2 (p = 1/3) for the collapsing phase of PBB. 
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IV. INCLUSION OF HIGHER-ORDER CORRECTIONS IN THE STRING FRAME 

-DILATON-DRIVEN CASE 



In the context of PBB cosmology, the natural frame to use in order to define the correction term Cc in the Lagrangian 
is the string frame. In this case, we should use / = e~'^R and u) = —e~'^, i.e. 



-0 



1 



+ JO.C 



(4.1) 



where (J) corresponds to the dilaton. This Lagrangian was suggested in [26,27], and used to construct nonsingular 
background cosmological solutions of PBB cosmology in the absence of a potential for the dilaton. Fluctuations in 
this model were studied in [36] in the absence of a dilaton potential. 

In the following, we extend these analyses to the case of a non-vanishing dilaton potential. We will first construct 
nonsingular solutions which in the Einstein frame correspond to nonsingular bouncing cosmologies. We then study 
how the fluctuations evolve across the bounce and compute the spectrum of fluctuations. The analysis in this section 
thus applies immediately to the modified PBB scenario in which the dilaton has a negative exponential potential. 
We can also apply the results to the initial version [5] of the Ekpyrotic scenario. In this case, the brane collision 
occurs at if = 0, and since thus the gravitational coupling constant does not change significantly near the bounce, the 
difference in the role of the higher derivative terms between the Einstein frame (the frame in which it appears most 
logical to define the correction terms in the Lagrangian) and the "string frame" (quotation marks used here because 
in the case of the Ekpyrotic scenario is not the dilaton, the dilaton being fixed) is not expected to be significant. 
The application of the results of this section to the version of the Ekpyrotic scenario with moving boundary brancs 
[38] and to the cyclic scenario [39] is more problematic since in this case ip is the dilaton, ip — oo at the bounce, 
and thus the difference in the evolution of models with correction terms defined in the Einstein and string frames is 
expected to be important. 

The reason why nonsingular solutions are possible in the presence of the correction term Cc, is that such a term 
can lead to violations of the null energy condition (from the perspective of an observer using unmodified Einstein 
equations). Thus, it is expected to lead to a successful graceful exit, in the same way that introducing matter violating 
the null energy condition allowed the construction of nonsingular bouncing models in [46,47]. 

In this model the background equations are written as 



6H'^ - 6H(j) + 4>^- 2Vs = efpo , 
A(j)H -AH- QH'^ -^"^ + 24) + 2Vs = e'^Pc , 
6H + 12ff2 + ^2 _ 2^ _ 6iJ^ _ 2{Vs - V^) = e't'A^ . 



(4.2) 
(4.3) 
(4.4) 



The dilatonic corrections Cc are the sum of the tree-level a' corrections and the quantum n-loop corrections (n 
1, 2, 3, • • •), with the function ^{^) [see (3.2)] given by 



C(.^) = -^C„e("- 



(4.5) 



71=0 



whore C„, (n > 1) arc the coefficients of ri-loop corrections with Co = 1. In this case the source terms due to Cc on 
the right hand side of (4.2)-(4.4) are given by [42] 



(4.6) 



71=0 



71=0 



71=0 



where 



{pc}^ = a'Af/^e^"-!)-^ {-24c(n - 1)H'^ + S^^} , (4.7) 
{Pc}„ = a'Ae("-i)* {8c(n - 1)H [(n - 1)^''H + 4>H + 2^{H + H^)^ + ^ (4.3) 
{A^}„ = a'Ae^"-^)"^ {24c(n - 1)H'^{H + H'^) - 3#^ [40 + A^H +{n- 1)0^] } , (4.9) 
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with A = —1/4. Following ref. [36] wc choose the coefficients as c = —d = —1. Note that the above corrections include 
the a' corrections (3.6)-(3.8), corresponding to n = 0. 

It is also convenient to relate the Hubble parameter and its derivative in the Einstein frame with those in the string 
frame by using eq. (2.10) 

= e*/^ (hs -fj, H^ = e^(Hs-^ + ^^Hs - ^ j . (4.10) 

Here the dots on the right hand side denote the time-derivative with respect to ts- The energy density pE and the 
pressure pe in the Einstein frame are expressed as 

PE^'iHl, pe---SHI-2He. (4.11) 

Once we know the evolution of the background in the string frame, it is easy to find the evolution of He, aE,f = —(j)/V2 
and to check whether the null energy condition, ps + Pe > 0, holds or not in the Einstein frame by using eqs. (4.10) 
and (4.11). Note that in the absence of higher-order corrections (£c = 0) one has 2He = —{pe +Pe) = —4>^ < 0. In 
this case once the contraction begins {He < 0) the Hubble parameter is always negative. Therefore it is not possible 
to have bouncing solutions required for the nonsingular Ekpyrotic scenario unless higher-order corrections Cc are 
taken into account. 



A. Background evolution 

In the absence of a negative exponential potential (Vg = 0), it was found in Ref. [26] that curvature singularities can 
be avoided by taking into account higher-order corrections Cc- In this case we have nonsingular bouncing solutions in 
the Einstein frame due to the violation of the null energy condition. We are interested in whether singularity avoidance 
is possible or not in the presence of the Ekpyrotic potential (2.9). Note that since near the bounce H'^ higher 
curvature corrections to the Einstein action will likewise be important in the presence of a potential. 

When Vs and Cc = the background solutions arc described by eqs. (2.8) and (2.12). In the string frame 
the scale factor evolves super- inflationary with growing Hubble rate (Hs > 0). We plot in Fig. 1 the evolution 
of background quantities both in the string and Einstein frames [see the case (i)]. The dilaton (p starts out from 
the weakly coupled regime g string = e"^ <C 1, corresponding to widely separated branes in the Ekpyrotic scenario, 
(p = — (/)/-\/2 ^1. In the Einstein frame the Universe is contracting with a negative Hubble rate. The solution 
inevitably meets a curvature singularity as (j) grows toward the strongly coupled regime (g'^tj-ing ^ !)• 

Our first main finding is that with Vs{4>) ^ there exist nonsingular trajectories in the presence of higher-order 
corrections {Cc ^ 0). Thus, the presence of the potential for the dilaton does not prevent the higher derivative terms 
from being able to smooth out the curvature singularity. The details depend on the value of the power-law index, p. 
When p ^ 1 the Ekpyrotic potential (2.6) is exponentially suppressed for </? > 1, in which case the dynamics of the 
system is hardly affected by the negative potential except for the region, ~ 0. However, in this region the higher 
derivative terms play a crucial role. 

In our simulations, we have adopted the potential (2.6) for > and Ve = for ip < 0. This is in the spirit of 
the first version of the Ekpyrotic scenario [5] in which the potential vanishes at the brane collision, the bulk brane is 
absorbed by the orbifold fixed plane via a small instanton transition, and there is no potential left afterwards. We 
show in Fig. 1 the dynamical evolution of the system for p = 0.1. The case (i) is the one in which only tree-level 
terms are present and in which singularity avoidance is not possible. The case (ii) corresponds to the one where both 
tree- level and one- loop corrections are taken into account (Ci = 1.0 and C2 =0). Inclusion of one-loop corrections 
makes it possible to have nonsingular cosmological solutions. In fact pE + Pe becomes negative around ts ^ 115 in 
Fig. 1, after which the Hubble parameter He begins to grow. The Universe starts to expand once He crosses zero. 
Namely the violation of the null energy condition allows to have nonsingular bouncing solutions in the Einstein frame. 
Nevertheless we should notice that the scale factors evolve super-inflationary both in the string and Einstein frames 
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due to unbounded increase of Hs and He, together with rapid growth of the field (t>. Therefore we are faced with 
another problem, namely how to connect to the stage of a decreasing Hubble parameter. 

If two loop terms are added (keeping the previous tree-level and one loop terms) phenomenologically more appealing 
nonsingular solutions can be obtained. When C2 is positive, the evolution of the system does not differ significantly 
compared to the case (ii). However, it is possible to obtain a decreasing Hubble rate if we take a negative value of C2. 
The case (iii) of Fig. 1 corresponds to the coefficients Ci = 1.0 and C2 = —1.0 x 10^''. We find that the growth rates 
of the scale factor and of are slowed compared to the case (ii). We see that ps + Pe becomes positive and begins 
to decrease toward +0 after the short period of violations of the null energy condition. Although this case does not 
correspond to the radiation-dominated Universe after the graceful exit, it is possible to connect to it by taking into 
account the decay of the dilaton to radiation ^. However, including radiation in the Ekpyrotic cosmology has some 
subtle points, and we do not consider this problem in the present work. 

We have checked that the addition of three loop terms with coefficients chosen to be of the order 10~^ (roughly the 
same hierarchy of coefficients between the two and three loop terms as between the one and two loop terms) does not 
change the results of the two loop analysis in a significant way. With a coefficient of the three loop term of the order 
1, the backgroimd solution ceases to be nonsingular. 

We emphasize that we have nonsingular bouncing solutions in the Einstein frame even in the presence of a negative 
exponential potential. When p <^ 1 the potential is vanishingly small fov ip ^ 1, in which case the dynamics of the 
system is practically the same as that of the zero potential discussed in Ref. [27]. In this case the dilaton starts out 
from the low-curvature regime \(p\ ^ 1, which is followed by the string phase with linearly growing dilaton and nearly 
constant Hubble parameter. During the string phase one has [26] 

(if 

asoc(-?7s)"\ (/> = --/-ln(-r?s) -l-const, (4.12) 

where </>/ ~ 1.40 and Hf ~ 0.62. In the Einstein frame this corresponds to a contracting Universe with 

aE oc {-VE)^f/^^"f^-' . (4.13) 

On the other hand, we can consider the scenario where the negative Ekpyrotic potential dominates initially but 
the higher-order correction becomes important when two branes approach sufficiently. Numerically we confirmed that 
it is possible to have nonsingular solutions (see Fig. 2). In the simulations we included the correction terms of £c 
only for (p < 1. In this case the background solutions arc described by eq. (2.8) or (2.12) before the higher-order 
correction terms begin to work. Given this background solution, one can obtain the spectra of curvature perturbations 
analytically, as we will see in the next section. The spectra depend on whether the higher-order terms are always 
dominant or not relative to the negative potential before the bounce. 



If wo were to include production of radiation at a fixed time during the expanding pliase, we could use the well-known results 
on the constancy of TZ in the expanding phase [44,48,49] to argue that the spectrum of fluctuations on large scales will be the 
same as what is obtained in this paper. The crucial fact about our bounce is that it is not symmetric in time (see Fig. 1). 
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FIG. 1: The evolution of Hs, He, as, as, 4>, and ps + Pe with c = —1, d = 1, p = 0.1. 
We choose initial conditions = —20, H = 5.0 x 10~^. The cases correspond to (i) only tree- 
level correction terms but no higher-order corrections (Ci = C2 = 0), (ii) tree-level and one-loop 
corrections present (Ci = 1.0, C2 = 0) (iii) tree-level and one- and two-loop corrections present with 
Ci = 1.0 and C2 = -1.0 x 10"^. 
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FIG. 2: The evolution of ip and as with c = —1, d=\,p = 0.1, 
Ci = 1.0 and Ci = -1.0 x 10"^. In this CEise we include the 
correction term Cc only for < 1. We choose initial conditions 
(j) = —15, H = 1.5 X 10~^. Prior to the collision of brancs at 
ip = 0, the universe is slowly contracting, which is followed by 
the bouncing solution through higher-order corrections. 



We have also studied Ekpyrotic potentials with other values of p, and found that if the potential is negative, 
corresponding to p < 1/3, then singularity avoidance is possible for suitable choices of Ci and C2 as in the case 
p = 0.1 shown in Fig. 1. In the case of p > 1/3 the field (p climbs up a positive exponential potential due to the 
Hubble contraction term. When p > 1/2 with Vq = p{l — 3p) we found that the field (f) returns back before it reaches 
the strongly coupled region, (f> > 0. This is equivalent to the fact that two parallel brancs do not approach each 
other sufficiently. In such cases the positive exponential potential makes the field bounce back before the higher-order 
correction becomes important (this may be related to the instability discussed recently in [50]). If we choose smaller 
values of Vb, it is possible to have nonsingular bouncing solutions which are similar to Fig. 1. This case corresponds 
to the one where the effect of the positive potential is negligible compared to higher-order corrections, in which case 
the backgromid solutions are given by cq. (4.12). When the positive potential is dominant from the beginning, it is 
difficult to obtain a solution where a successful graceful exit is realized by higher-order corrections. 
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B. Density perturbations 



Let us proceed to the analysis of the evolution and the spectra of density perturbations. We shall consider two 
cases: (i) the effect of the potential V{(p) is always negligible relative to the correction term Cc, and (ii) the effect of 
the correction term Cc becomes important only around the graceful exit (0^0). Note that the second case is the 
physically more interesting one for applications to Ekpyrotic cosmology. 

1. Case (i).- \V{4,)\ < \£c\ 

When the correction terms (3.6)-(3.8) always dominate relative to the exponential potential (2.6), the spectra of 
density perturbations are similar to the ones discussed in ref. [36]. During the string phase with linearly growing 
dilaton and nearly constant Hubble parameter with (j)f ~ 1.40 and Hf ~ 0.62, we have a sufficient amount of inflation 
with e-folds N = \n{a/ai) > 60 provided that the dilaton field satisfies |</)| ^ 1 initially [36,42]. In this stage Q defined 
in eq. (3.13) is proportional to e~'^ by making use of eq. (4.12), thereby leading to 

z(x{-r]sy, with 7 = -1 + ^-0.13. (4.14) 

Making use of the relation (3.21) which is valid for positive s, the spectral tilt of the large-scale curvature perturbation 
is 

~2.26. (4.15) 

The evolution of the frequency shift s is nontrivial (see Fig. 3). In the low curvature regime where the higher-order 
terms are not important, s is positive (s ~ 1), as in the usual PBB scenario. It then changes sign and becomes 
negative during a short transition from the low-curvature regime to the string phase. During the string phase, ^ and 

Hs arc constant (0 ~ 1.40 and Hs — 0.62), and ^ ^ —e~'^. The correction term on the right hand side of eq. (3.13) 
for s dominates in this phase. It follows from eq. (3.13) that the (/) dependence of the leading term cancels out between 
and u){ip), and that hence s is constant and negative until the graceful exit. In a stage with negative constant s, 
the solution of eq. (3.14) can be written in the form 

= VW\ [ci/.(x) + C2K,{x)] , (4.16) 

where x and v are given in (3.17) (with s replaced by the absolute value of s), and /,y and are modified Bessel 
functions, whose asymptotic solutions are I^, oc x" , oc x~'' for x ^ 0, and ~ /^/2'kx, K^, ~ ^7r/(2a;)e~^ 
for a; — > oo. Then one reproduces the spectral tilt (4.15) in the large-scale limit (|sfc^| <C \z"/z\). For small-scale 
modes curvature perturbations show exponential instability due to negative frequency shift. After the horizon crossing 
([sfc^l < \z" /z\), curvature perturbations are frozen, since 7 is smaller than 1/2 in this case. 

It was shown in ref. [31] that the ratio (j)f/Hf is required to lie in the range 2 < (f)f/Hf < 3 for the successful 
graceful exit in the presence of other forms of higher-order a' correction. Therefore the spectral tilt lies in the range 

2<nTC-l<3, (4.17) 

which is valid for large-scale modes (|sfc^| <C |2:"/^;|). Therefore we have blue-tilted spectra as long as the correction 
jCc dominates compared to the exponential potential. 
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FIG. 3: The evolution of the frequency shift s for c=— l,d=l,p = 
0.1 with initial conditions = —100, H = 1.5 x 10~^. We include the 
quantum correction £c from the beginning. The shift s is approximately 
constant (and negative) during the string phase, which is followed by 
the stage of decreasing curvature with positive s. Inset: The evolution 
of s in the case where the quantum correction is taken into account 
only for (/? < 1. Note that s rapidly changes sign around the graceful 
exit. 



2. Case (ii) ; \V{4>)\ > \La\ but < \Lc\ for ~ 

When the correction term £c becomes important only around the graceful exit (0 ~ 0), the spectra of density 
perturbations generated before the bounce are mainly determined by the exponential potential. In this case the 
evolution of the background can be characterized by eq. (2.12). Then the quantity Q in eq. (3.12) evolves as 

^ 20^e-^ ^ (_ )2VP/(i-P) . (4.18) 
(2i/s-<^)2 

Therefore we find 

zcx(-r,s)\ with 7=7^, (4.19) 

1 -p 

and the spectral tilt for the curvature perturbation is 

/ (for < p < 1/3) , 

(lorl/3<p<l). ^'-^'^ 

This coincides with the result in the Einstein frame performed in ref. [14-17]. For very slow contraction with a 
negative Ekpyrotic potential (p <C 1), one has blue tilted spectra with titi — 1 = 2. Since 7 is less than 1/2 for p < 1/3 
(i.e., negative potential), curvature perturbations are not enhanced in the large-scale limit even in the presence of the 
correction Lc around the graceful exit. The simplest PBB scenario with zero potential corresponds to p = 1/3 and 
7 = 1/2, in which case one has n-ji — 1 = 3. In this case curvature perturbations evolve as TZk oc ln(— 77) as found from 
eq. (3.15) with eq. (3.16). 

We have solved the evolution equation (3.11) for the cosmological fluctuations numerically. Experience from study- 
ing fluctuations in inflationary cosmology teaches us that following the evolution equation for 5" instead of for the 
gravitational potential $ is less likely to be effected by numerical noise. Since in a contracting Universe, one of the 
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modes of $ increases much more rapidly than the dominant mode of '5. we beheve that it is advantageous to use ^! 
in our case as welL In addition, from a more conceptual point of view, the variable TZ is preferable since it is more 
closely related to ^ in terms of which the action for cosmological fluctuations takes on its canonical form. Note that 
^ is also the good variable to use when following cosmological fluctuation from inflation through reheating [70]. In 
Fig. 4 we plot the resulting evolution of the spectra of curvature perturbations for several different frequencies. The 
higher-order correction Cc is included only when two branes approach sufficiently, i.e., ip < I. We find that large-scale 
modes {k <C 1) are not enhanced as predicted by eq. (3.15). In contrast, small-scale curvature perturbations exhibit 
rapid growth around the graceful exit 0). 

There are two reasons for this instability. The first is the fact that the frequency shift s becomes negative for a 
short period where the higher-curvature effect is dominant (see the inset of Fig. 3). As is obvious from (3.14), an 
exponential instability for "if is induced by negative s, which is stronger for larger k. We expect this instability will 
become important for modes with |s|i„ax k^ > 1, where |s|inax is the maximal negative value of the function s. Prom 
the inset of Fig. 3 the maximal absolute value of s during the negative branch is about |s|max ^ 10^. Hence, the 
s-instability is expected to be important only for modes with k > 10~^. By comparing runs with s given by the 
general formula and runs with s = 1, we were able to determine numerically that the actual cutoff value of k below 
which the instability due to the s-term is negligible is fc ~ 10"^. Thus, we conclude that the main source of the 
short wavelength instability of the fluctuation modes around the bounce must be a second one, namely the nontrivial 
nature of the bouncing background and its result on the quantity z"/z. 

After the transition to the expanding Universe, the curvature perturbation is nearly conserved as found in Fig. 4. 
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FIG. 4: The evolution of the spectra of curvature perturbations, Pjz, for 
c = —1, d= l,p = 0.1, Ci = 1.0 and C2 = -1-0 x 10"^. The initial conditions 
are chosen to be = —15, H = 1.5 x 10~^. We include the higher-order 
correction Cc only ior ip < 1. The curvature perturbation does not exhibit 
significant variation during the contracting phase. However small-scale modes 
are enhanced around the graceful exit. 

We show in Fig. 5 the spectra of curvature perturbations for p — 10^'^ in the case where the correction £c is included 
only for (fi < 1. We find that the numerical value of the spectral tilt is — 1 2 for fc < 10~*, which coincides with 
the analytic estimation (4.20). However this estimate is no longer valid for small-scale modes due to the negative 
frequency shift and the instability around the graceful exit. The spectra are highly blue-tilted for k > 10^^ as found 
in Fig. 5. This growth of small-scale fluctuations obviously works as the gravitational back-reaction to the background 
evolution. Although we did not consider the effect of the back-reaction here, it is certainly of interest to investigate 
how the background evolution is modifled around the graceful exit. We have performed the simulations with various 
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choices of time steps to make sure that the effects we find are not numerical artifacts. The spectra obtained are 
independent of the specific value of the time step. 




FIG. 5: The final spectra of the curvature perturbation, P-jt, and of 
the gravitational potential, P#, in a simulation with c = —1, d = 1, p = 
10~^, Ci = 1.0 and C2 ~ —1.0 x 10~^ with same initial conditions as 
in Fig. 4. The superscript "b" denotes the quantities before the bounce 
(at ts = 300), while "a" indicates the quantities after the bounce (at 
ts = 1000). Wo have included the quantum correction £c for values 
1^ < 1. The spectral tilt is n-jz — 1 2 for A: < 10"'', which agrees 
with the analytic estimation of eq. (4.20). For the modes k > 10~*, the 
spectra arc highly bluc-tilted, due to an instability of small-scale modes 
during the graceful exit. The fluctuations in $ are nearly scale-invariant 
on large scales before the bounce. 



In Fig. 5 we have also plotted the induced fluctuations of in the Einstein frame, determined from the results 
for TZ^ in the string frame, using the relation 

^^<xnl/k\ (4.21) 

This corresponds to eq. (6.16) in Appendix B, which follows from the relation (6.14) in the Einstein frame. Note that 

this relation is valid in the absence of higher curvature corrections to the Lagrangian, and will therefore be good at 
times long before and long after the bounce. For the negative Ekpyrotic potential (0 < p < 1/3), one has < 7 < 1/2 
from eq. (4.19), in which case the second term in eq. (6.17) completely vanishes. In this case we have the relation 
(6.19), namely 

oc Hi'jf/k. (4.22) 

Note that H^]^^^ can be written by the sum of two terms which arc proportional to {k\ris\Y~^ and {k\r]s\)~'^~^^- 

Since < < 1/2 for < 7 < 1/2 (i.e., < p < 1/3), the term proportional to {k\r]s\y~^ is the growing mode 
during the contracting phase on large scales. Therefore the spectrum of $f before the bounce can be estimated as 

P| DC fc^'^-i DC fc-^T oc , (4.23) 

from which we have 

n$-l = ---^. (4.24) 
I— p 
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Then \vc hiwc a scale-invariant spectrum before the bounce for p ^ 0, as first pointed out in [18]. This agrees with 
our numerical result shown in Fig. 5. 

The term proportional to {k\r]s\Y^^, however, decays after the graceful exit as long as < 1. The dominant mode 
in $1' long after the bounce is described by the term {k\ri\)~'^~^^ , in which case the spectrum of $ is written as 

P| cx k^-^" . (4.25) 

From Fig. 5 wc find that the spectrum of <I> is blue-tilted with n<j, ^ 3 for k > 10^^° (small-scale modes for k > 10~^ 
exhibit larger blue tilt with n$ > 3). This corresponds to the value i/ ~ 1/2 in eq. (4.25) after the bounce. 

Our numerical calculations show that large-scale modes with k < 10~^° do not exhibit such a blue spectrum. 
This can be understood that the term proportional to {k\r]s\Y~^ which is dominant in the contracting phase does 
not become smaller than the one proportional to (fc|?7s|)~''~''^ in the expanding branch for very small k, unless we 
evolve the fluctuations until long after the bounce. However it is rather difficult to follow such large amount of time 
numerically. In addition the second term in eq. (6.17) is not numerically negligible relative to the first term for these 
large-scale modes due to the modification of the equation of state after the bounce [when 7 > 1/2 the second term in 
eq. (6.17) is nonvanishing as found by eq. (6.18)]. Nevertheless, we expect that the term proportional to {k\'r]s\)~'^~^^ 
in the first term in eq. (6.17) eventually dominates long after the bounce, in which case the spectrum is given by 
eq. (4.25). Therefore the final spectrum of $ is not generally scale-invariant. The spectral index is dependent on 
the evolution of the background after the bounce (i.e., 7). In this sense including radiation is necessary in order to 
evaluate the spectrum of $ in realistic cases where the solution connects to our Friedmann branch. 

From Fig. 6 we find that the amplitude of $ decreases after the bounce, thus showing that the dominant prc-bouncc 
mode of $ couples exclusively to the decaying mode of $ after the bounce, as derived in [15] using matching conditions 
on a constant scalar field hypersurface. 
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FIG. 6: The evolution of the curvature perturbation, Pr., and of the 
gravitational potential, Pg, for the fluctuation mode corresponding to 
k = 10~^. The model parameters and initial conditions are the same as 
in Fig. 5. The amplitude of the gravitational potential near the bounce 
when the higher derivative terms dominate cannot be trusted since $f 
is computed from TZf via the equation (4.21) which is only valid in 
the absence of such higher derivative terms. As follows from this plot, 
the dominant growing mode of $f during the period of contraction 
couples only to the post-bounce decaying mode. At the time of the 
bounce, the curve for $ is dominated by numerical noise. However, 
since $ is computed at each time separately from the value of TZ at 
that time, this does not introduce numerical errors in the late time 
values of 3>. 
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Eqs. (3.20) and (4.19) indicate that a scale-invariant spectrum may be obtained for p = 2/3 for the modes which are 
enhanced during the bouncing phase [23]. In order to obtain such a spectrum, the exponential potential (positive in 
this case) is required to dominate the higher-order term except around the graceful exit. However we have found that, 
for some likely initial conditions, the field bounces back toward larger before the higher-order correction begins to 
work. In language appropriate to Ekpyrotic cosmology, this means that the branes never collide. If the higher-order 
term always dominates compared to the positive exponential potential, we have blue-tilted spectra (4.15). 

V. INCLUSION OF HIGHER-ORDER CORRECTIONS IN THE EINSTEIN FRAME 

-MODULUS DRIVEN CASE 

In this section we consider adding higher derivative terms defined in the Einstein frame. We add a Gauss-Bonnet 
term proportional to i?QB multiplied by a function of the modulus field ip to the action. Such a term arises as the 
one loop correction in the context of orbifold compactifications of the heterotic superstring [51]. Since the initial 
version of Ekpyrotic cosmology [5] is based on an orbifold compactification of a theory dual to heterotic superstring 
theory, the correction terms used in this section are well motivated in the context of the scenario of [5]. Indeed, it 
was found in Rcf. [52] (see also [53]- [60]) that the inclusion of the Gauss-Bonnet term coupled to a modulus field in 
the Einstein frame leads to the possibility of obtaining nonsingular solutions. In the work of [52], the potential for 
the modulus field was taken to vanish. In this section we will include an exponential potential.^ More specifically, 
the correction Lagrangian we consider here corresponds to [in the notation of (3.1) and (3.2)] f = R, oj = 1, c = — 1, 
d = 0, ^{(fi) = ln[2e'^r]'^{ie'^)] with r]{ie'^) being the Dedekind ry-function [52]. Here ^{(p) is approximately given by 

^(<^)^-^(e^ + e-^) . (5.1) 

The sign of A is chosen to be positive, which is different from the one discussed in the previous section. Note also 
that even though the coefficient ^(<^) becomes large at large brane separation (large negative values of the dilaton in 
the case of PBB cosmology), this increase is outweighed by the falloff of the curvature invariant, as in the case of the 
model considered in the previous section. Thus, in Ekpyrotic cosmology the correction terms are expected to become 
important only in the high curvature region. 

Let us analyze the one-field system of a modulus tp, keeping the dilaton fixed. We will consider solutions starting 
in an asymptotically fiat region and beginning in the expanding branch. We have not found solutions which begin in 
a contracting phase and undergo a successful bounce. However, note that in the original Ekpyrotic scenario of [5], the 
scale factor on the orbifold fixed plane corresponding to our four-dimensional space-time corresponds to an initially 
asymptotically flat region, and is always expanding. Thus, the solutions found here might be applicable to a version 
of Ekpyrotic cosmology formulated entirely in terms of physics on the orbifold fixed plane. 

A. Background evolution 

As was discussed in [52], when Ve = the PBB singularity can be avoided for positive values of A when the 
a' corrections introduced above are taken into account. The sign of A is crucial for the existence of nonsingular 
cosmological solutions. For negative values of A, the a' corrections do not help to lead to a successful graceful exit, 
as was analyzed in ref. [59]. 

In the absence of the Ekpyrotic potential, the background evolution for ts < is given by [56] 



®The authors of ref. [60] analyzed nonsingular cosmological solutions in the presence of some positive potentials (not the 
exponential potential) . 
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0.E — ao, 



- -2-, 

''E 



(5.2) 



where ao and Hq (> 0) are constants. The Gauss-Bonnet term leads to a violation of the null energy condition 
ifE + Pe < 0) at sufficiently large curvatures and thus enables a graceful exit. If we start in an expanding branch 
(contrary to the spirit of PBB and Ekpyrotic cosmology), this leads to a super- inflationary solution {He > 0) until a 
"graceful exit" (see Fig. 7). The Universe is initially expanding very slowly with a nearly constant scale factor. After 
the Hubble parameter reaches its peak value He — H^i^-^, the system connects to a Friedmann-like Universe with 

1 /3 

He — l/{3tE), ciE , and (/J c!c — Ini^;. 

If the Ekpyrotic potential is present, the situation is quite different. We have adopted the potential (2.6) for ip > 
and Ve ~ for ip < 0. Once again, we start in an expanding phase. Initially, the potential term is not important 
and the Universe evolves in a super-inflationary trajectory until a graceful exit after which the Hubble expansion rate 
begins to decrease. When p < 1/3, corresponding to the case of a negative exponential potential, then as i/j ^ 0, 
the potential becomes important and leads to a change in sign of He- We find that the system enters a stage of 
slow contraction [see the case (ii) of Fig. 7]. Note that in Fig. 7 He changes the sign twice. After the negative 
Hubble peak, the Hubble rate begins to grow toward He ^ — without changing sign. Then the system enters a 
very slow contracting phase with a nearly constant scale factor. In this stage the field ip evolves rapidly toward large 
negative values. In the presence of negative exponential potential {p < 1/3) we have found that the contracting stage 
eventually appears even when Vq is small. 
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FIG. 7: The evolution of He, cle, ip, and pE + Pe in the modulus-driven case with c = —1, d = 0. We 
choose initial conditions <^ = 20 and H = 5.087 x lO"**. Each case corresponds to (i) Vo = with p = 0.1, 
(ii) Vq = 0.01p(l - 3p) with p = 0.1, (iii) Vo = p{l - 3p) with p = 0.5. 



When p > 1/3 there exists a positive potential barrier as the field (p approaches zero. The case (iii) in Fig. 7 
corresponds to p = 1/2 with Vq = p{l — 3p). The effect of the positive potential is important around ip ^ 0, which 
works to return the field back toward larger p. After the graceful exit the Hubble rate is always positive with slowly 
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changing (p. The scale factor evolves as a power-law (a cx fP) due to a positive exponential potential. The p > 1/3 
case provides us with reasonable nonsingular cosmological solutions. Nevertheless we need to caution that these 
nonsingular solutions are different from the bouncing ones where the contraction of the universe occurs before the 
graceful exit. 

One may argue that the bouncing trajectories may be found by including the correction Cc only around </? ~ 0. 
However we have numerically found that this is not the case. The super-inflationary evolution characterized by 
eq. (5.2) is typically required for the construction of nonsingular solutions in the present scenario. 

B. Density perturbations 

When the Gauss-Bonnet term is dominant relative to the Ekpyrotic potential, the spectra of density perturbations 
can be analyzed as in the case of the zero potential (p = or p = 1/3). In this case the evolution of the background 
during the phase of modulus-driven inflation is given by eq. (5.2), thereby leading to — —{TT/3)ipe'^ cx (— ^b)^. 
Making use of this relation together with eq. (3.12), we find the evolution of Q and z as 

Qcx(-tij)2, z cc {-tE) cc {-r]E) ■ (5.3) 

This means that 7 = 1 in eq. (3.16), in which case curvature perturbations are enhanced on super-Hubble scales 
during super-inflation {TZk oc (— ??e)~^) due to the growth of the second term in eq. (3.15) [61]. We show in Fig. 8 the 
evolution of curvature perturbations in the case of zero potential (p = 1/3) for two different modes {k = 10~^ and 
k = 10~^). We find that curvature perturbations are amplified before the graceful exit. 

In order to obtain the spectral tilt of density perturbations, we have to caution that the function s defined by 
eq. (3.13) is a time-varying function and is proportional to {—Ie)- Therefore the formula (3.21) can not be directly 
applied. Instead one is required to consider the evolution equation for curvature perturbations: 

2 ■ k'^ 

nu + —Uk - a—tEUk = , (5.4) 

lE Oq 

where a (> 0) is a constant that depends on Hq in eq. (5.2). The solution of this equation is written in terms of the 
Bessel functions 

Tlk = {-tE)-^'^ [ci J_i/3(a;) + C2 Ji/3(a;)] , (5.5) 

where x = |-\/a^(— i^)^^^. Notice that this solution asymptotically approaches the Minkowski vacuum for a; — > 00. 
Since J±i/3(a;) a fc^^/^ in the a; — > limit, the spectrum of large-scale curvature perturbation is proportional to 
P-jz oc k"^/^. Therefore the spectral index is 

7 

nn-l = (5.6) 

which is a blue-tilted spectrum. 

In the absence of the Ekpyrotic potential, the evolution of the background in the asymptotic future is given by 

I 1/3 1/3 \ / '2, 

ip ~ -\/3/21ni£;, He cx 1/(3^^;), and acct^ . Therefore one has z oct^ oc t]^ in eq. (3.15), in which case curvature 
perturbations exhibit logarithmic growth, 

TZ (X In rjE- (5.7) 

This indicates that the second term in eq. (3.15), which we call "D-mode", dominates even after the graceful exit. In 
the case where the D- mode decays after the graceful exit, the surviving spectra observed in an expanding Universe 
should correspond to the first term in eq. (3.15) ["C-mode"]. In the present model, however, the D-mode survives in 
an expanding branch. Therefore the spectrum of the curvature perturbation during super-inflation can be preserved 
even after the graceful exit. In fact, the numerical value of the final spectral tilts of TZ are found to be n-jz — 1 2.3 
for the modes fc <C 1 (see Fig. 9). This agrees well with the analytic result (5.6). 
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FIG. 8: The evolution of the spectra of curvature perturbations 
in the modulus-driven case (c = —1 and d = 0) for {p,k) = 
(1/3, 10"^), (1/3, 10"=*), (1/2, 10"^), (1/2, 10"^). Wo choose initial con- 
ditions ifi = 23.888 and H = 4.158 x 10"". Note that the p = 1/2 
case corresponds to the positive exponential potential while the zero- 
potential corresponds to p = 1/3. Around the graceful exit curvature 
perturbations exhibit rapid growth especially when the potential is pos- 
itive (p > 1/3). 




FIG. 9: The final spectra of curvature perturbations in the modulus- 
driven case (c = —1 and d = 0) for p = 1/3 and p = 1/2. The initial 
conditions are the same as in Fig. 8. When the positive potential is 
present (p > 1 /3) , the amplitude of the spectrum is larger than in the 
case of the zero-potential. We numerically find that the spectral tilt is 
n-R. — 1 ~ 2.3 for fe <C 1, which agrees well with the analytic estimation, 
mi-l = 7/3. 



When a positive Ekpyrotic potential is present {p > 1/3), the dynamics of the system is more unstable around 
the graceful exit. This leads to the violent growth of curvature perturbations when the field bounces back due to 
the potential barrier. This threatens the viability of the cosmological perturbation theory around the graceful exit. 
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Nevertheless the perturbations are not singular as long as the background is smoothly joined to the expanding branch. 
We have numerically evaluated the power spectra of TZ for the modes which left the horizon during super-inflation. 
Although the amplitude is larger compared to the case of the zero-potential (p = 1/3), the final spectral tilts are 
similar for large-scale modes {k <^ 1), see Fig. 9. Again the final spectra are found to be blue-tilted. We should also 
mention that the frequency shift s becomes negative for the Hubble rate which is larger than unity around the graceful 
exit [61]. In this case the small-scale modes show exponential instability as we pointed out in the dilaton-driven case. 
The negative Ekpyrotic potential (p < 1/3) is not worth studying, since this case does not connect to the expanding 
branch as analyzed in the previous subsection. 

Finally, we should mention that we have neglected the effect of radiation in all our analysis. However this is 
expected to appear at some moment of time. This can also alter the final spectra of curvature perturbations due to 
the dominance of the C-mode in eq. (3.15). We leave to future work about these realistic situations. 

VI. DISCUSSIONS AND OPEN ISSUES 

We have studied the effects of higher derivative terms in the joint gravitational and matter action for theories 

motivated by Pre-Big-Bang and Ekpyrotic cosmology with a single scalar matter field with an exponential potential. 
Applied to PBB cosmology, our model corresponds to a theory with an exponential potential for the dilaton. In the 
language of the initial version of Ekpyrotic cosmology [5] , our scalar field is the modulus field corresponding to the 
separation of the bulk brane from our orbifold fixed plane, in the second version [38] and in its cyclic version [39,40] 
the field is the radius of the extra spatial dimension. The higher derivative terms introduced are the leading string 
and quantum corrections to the low energy effective action of string theory. 

When applying the correction terms in the string frame, and for suitable choices of the coefficients of the higher-order 
corrections we find nonsingular cosmological solutions which in the Einstein frame correspond to bouncing Universes. 
We thus find that higher derivative terms can smooth out the singularities in PBB and Ekpyrotic cosmology and 
lead to a graceful exit (in the language of PBB cosmology) or a nonsingular bounce (in the language of Ekpyrotic 
cosmology). We have thus generalized the results of [26] to models with exponential scalar field potentials. 

We have studied the evolution of fluctuations in our nonsingular bouncing cosmologies. This analysis is not plagued 
by the matching ambiguities inherent to analyses where the contracting and expanding cosmologies are matched across 
a singular space-time surface. For all potentials with < p -C 1 we find a blue spectrum of curvature fiuctuations. 
The precise spectral index depends, as expected, on whether the higher derivative correction terms arc important 
at times when the scales on which we compute the fluctuation spectrum exit the Hubble radius during the phase of 
contraction. 

If the higher derivative terms arc not dominant when the scales exit the Hubble radius, the index of the spectrum 
agrees with what is obtained by applying the general relativistic matching conditions on a uniform density hypersurface 
[15-17]. The only difference is an instability of small-scale fiuctuation modes during the bounce (see also [61]) which 
leads to a further steepening of the spectrum. Our result implies that the growing mode of $ during the contracting 
phase which is scale-invariant for < p <C 1 is effectively uncoupled with the dominant constant mode of $ in the 
expanding phase, a result obtained in the context of matching conditions in [13] (in the case of PBB cosmology) and 
in [15,16] for the Ekpyrotic scenario. If the higher derivative terms dominate when scales of interest exit the Hubble 
radius, then the spectrum is blue with a slope of — 1 ~ 2.26. Note that our result implies that it is the curvature 
fluctuation TZ (more precisely the variable C originally introduced by Bardeen in [62] and used in [15], which equals 
TZ up to terms which are suppressed by fc^ for large-scale fluctuations) which is effectively conserved for large-scale 
perturbations across the bounce. 

We have also studied nonsingular cosmological models obtained by adding a Gauss-Bonnet term (defined in the 
Einstein frame) multiplied by a suitably chosen function of the single scalar matter field in the problem (a modulus 
field). Once again, we have included an exponential potential for the modulus field. Although we do not find 
bouncing cosmologies, we find interesting nonsingular cosmological solutions which begin in an asymptotically flat 
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region, undergo super-exponential inflation followed by a graceful exit to a phase with decreasing Hubble radius. In 
the presence of a negative exponential potential (0 < p < 1/3), the solutions reach a maximal radius and begin to 
contract as the field crosses = 0. During this period of contraction, the Hubble parameter remains finite. Such 
solutions might be applicable to Ekpyrotic Universe models formulated in terms of physics on the four-dimensional 
orbifold fixed plane corresponding to our visible space-time. When the potential is positive {p > 1/3), the modulus 
(f bounces back around the brane collision toward larger ip due to the barrier of the positive potential. Although 
singularities can be avoided around ~ 0, this model does not correspond to bouncing solutions where the contraction 
of the Universe occurs due to the Ekpyrotic potential before the graceful exit. 

We have also studied the spectrum of curvature fluctuations TZ in this modulus- driven cosmology. When higher- 
order corrections are important before the bounce (as it must be for the existence of nonsingular solutions), one has 
n-jz — 1 ~ 7/3. This result is again in agreement with what can be obtained by neglecting the graceful exit and 
matching two Einstein Universes at a constant density hypersurface. 

Note that we have chosen to evolve the curvature perturbation TZ on comoving hypersurfaces, and found that 
the spectral index is given by (4.20) when the Ekpyrotic potential is dominant. Note that this spectrum in the 
< p < /1/3 case comes from the C-mode in eq. (3.15), which is blue-tilted for very slow contraction (0 < p <C 1). 
Since the large-scale D-modes are enhanced for 1/3 < p < 1 during the contracting phase, the spectrum of TZ will be 
scale-invariant for p 2/3 right after the bounce [23]. 

If we follow instead the gravitational potential $ in the longitudinal gauge, its spectral index generated during 
the collapsing phase is estimated as eq. (4.24), which is different from that of TZ [sec eq. (4.20)]. When p ^ 0, 
corresponding to a very slow contraction, the growing mode (D-mode) of $ is approximately scale-invariant. The 
authors of ref. [18,19] claimed that a scale-invariant spectrum of the dominant post-bounce mode of $ would inherit 
this scale-invariant spectrum. 

However, we know that when computed at late times long after the bounce, in an expanding universe, the spectra 
of TZ and $ must be identical. Thus, given our results concerning the spectrum of TZ, we know that the spectrum 
of $ long after the bounce cannot be scale-invariant. Our numerical simulations show that the contribution from 
the pre-bounce D-mode decays after the system enters the expanding branch, and thus shows that the pre-bounce 
growing mode of $ couples exclusively to the post-bounce decaying mode. This results in a blue-tilted spectrum of 
<f> when evaluated long after the boimce (see Fig. 5). For very large-scales with k < 10^^°, we need to solve the 
equation of fluctuations up to sufficient amount of time in order to find the complete decay of the D-mode relative to 
the C-mode. In addition the second term in eq. (6.17) is not numerically negligible for very small k when 7 is greater 
than 1/2. Nevertheless the term proportional to {k\ris\)~"~^^ in the first term in eq. (6.17) eventually dominates long 
after the boimce (77s 00), thereby yielding the spectrum (4.25). Therefore the spectrum of $ long after the bounce 
is not generally expected to be scale-invariant, whose spectral index depends on the evolution of the background in 
an expanding branch. 

Since near the bounce the magnitudes of the two modes of $ and TZ differ by such a large ratio, we must worry 
about the possibility of numerical errors. In particular, if one were to follow the evolution equation for it would be 
difficult to ensure that numerical noise does not lead to an artificial coupling between the pre-bounce growing mode 
and the post-bounce dominant (constant) mode. We have checked that our results do not seem to suffer from a similar 
problem by repeating the simulations with different values of the time step, At. We did not find any dependence of 
the results within the range of time steps we have chosen (10~^ < At < 10~^). 

Let us compare our findings to results which have already appeared in the literature. As mentioned repeatedly, 
our results concerning the spectrum of fiuctuations obtained in the classes of nonsingular bouncing Universe models 
considered in this paper agree with the results of [15-17] obtained when removing the higher derivative correction 
terms (thus going back to a singular background) and matching the fluctuations on a constant scalar field matching 
surface. The results imply that the growing mode of $ in the contracting phase does not source the post-bounce 
dominant mode of Our results thus indicate that the conjecture of [18,19], namely that the growing mode of $ in 
the contracting phase (which in Ekpyrotic cosmology has a scale-invariant spectrum) should generically determine the 
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amplitude and spectrum of the dominant mode of <I> in the post-boimcc phase, is not vahd. As emphasized in [24] and 
[19], in the case of a singular background the spectrum of fluctuations in the expanding phase depends sensitively on 
the details of the matching conditions used. Since we have only used one class of ways to smooth out the singularity, 
the sensitive dependence on the matching surface might not have been completely eliminated, but might find itself 
reflected in a sensitive dependence of the final spectrum on the specific form of the correction terms in the action. 
We leave the study of this issue to future work. 

Our work indicates that it is difficult to obtain a scale-invariant spectrum of curvature fluctuations for a single 
field PBB or Ekpyrotic cosmology. However, in the case of Ekpyrotic cosmology there is the intriguing fact that the 
growing mode of the gravitational potential $ during the phase of contraction has a scale-invariant spectrum. To 
obtain a scale-invariant spectrum of $ and thus also of the curvature fluctuation TZ at late times in the expanding 
phase, a suggestion [18,19] was to non- trivially connect the growing mode of $ during the contracting phase with 
the constant mode in the expanding phase. We have shown that this does not occur in the single field case with our 
choice of correction terms to the action (needed to obtain a nonsingular bounce). 

Note that there are examples where a large growth of $ during the phase of contraction persists after the bounce 
(see e.g. [23,47,37]). A criterion for when this occurs has been proposed recently in [37]. The condition is that the 
relative variation of TZ over a Hubble time scale should be appreciable, i.e. the following relation [37] 

(6.1) 

should hold close to or right at the bounce. We use eq. (3.15) and restrict to the case of an exponential potential, in 
which case one has 

7t _ Dk{l-P)v for0<p<l/3, 

HTZ pa^{Ck + £>fc(-77)i-27) | 1^ for 1/3 < p < 1 . ^ ' ' 

The marginal case p = l/2> (important for PBB and also for Ekpyrotic scenario, in which the potential disappears 
close to the bounce) should be treated separately and leads to: 

HTZ Ck + DuH-kri) ' ^ ' ' 

Surprisingly enough, both the results (6.2) and (6.3) indicate that $ could never match to TZ nontrivially. For 
p < 1/3, n/HTZ as ;/ 0. For 1/3 < p < 1, iZ/HTZ ^ (1 -p)/33 ~ 0(1) as ^ 0. This latter case is interpreted 
as a variation of TZ rather than a change induced by $. Interestingly enough, if one takes seriously the ratio (6.1), the 
singularity at the bounce (i.e., if ry = is reached or not) does not matter in the impossibility of matching $ to TZ. 

Recently several authors [47,37] considered models of a bouncing Universe (realized in [47] by introducing matter 
violating the weak energy condition and in [37] by making use of spatial curvature in the background metric) in which 
TZ grows dramatically across the bounce and there is a coupling between the growing mode of $ in the contracting 
phase and the dominant mode of $ in the expanding phase. In this case, it may be possible to obtain a scale-invariant 
spectrum, as already realized in [23]. 

Although we have concentrated on the density perturbation in the single field scenario, the situation can be changed 
by taking into account a second scalar field [23,63,64]. The system of multi-component scalar fields generally induces 
isocurvature perturbations, which can be the source of adiabatic perturbations. In such a case the relation (6.1) 
could be satisfied, since isocurvature perturbations act as source term for TZ in addition to (&. In fact the authors 
of ref. [63] considered a specific two-field system with a brane-modulus <^ and a dilaton When the dilaton has 
a negative exponential potential with a suppressed Ekpyrotic potential for ^p, the entropy "field" perturbation can 
be scale-invariant if the model parameters are fine-tuned [63]. It was also pointed out in ref. [23] that the quantum 
fluctuation of a light scalar field (with a non-canonical kinetic term as studied in [65]) such as axion may lead to 
the flat spectra of isocurvature perturbations. If the correlation between adiabatic and isocurvature perturbations is 
strong, adiabatic perturbations may be scale- invariant. 
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Wc wish to stress that our work is not conchisivc. In particular, in order to fully evaluate the final power spcx'tra. one 
should solve the equations of motion for fluctuations in a nonsingular bouncing model including radiation. Important 
issues which should be investigated further include: 

• The final power spectra of the curvature perturbation are found to be blue-tilted for the nonsingular Ekpyrotic 
models we have considered, which rely on specific higher derivative correction terms. Are there other correction 
terms to the action which are motivated by string theory, lead to nonsingular bouncing scenarios and yield a 
flat spectrum even in the single fleld case? Perhaps toy bouncing models using exotic scalar fields or matter in 
refs. [46,47] can be a good starting point to construct viable nonsingular Ekpyrotic models. 

• We do not include the effect of radiation (or particles) which can be efficiently produced near the bounce. 
Particle production around the transition region is expected to be quite efficient [66], and this could lead to an 
additional instability of small-scale metric perturbations. This effect may also non-trivially alter the nonsingular 
background evolution by the back-reaction effect of created particles. It is also required to include the radiation 
after the bounce in order to evaluate the surviving spectra accurately, although the coupling between the scalar 
field and radiation should be chosen carefully in that case. 

• It is of interest to study the effect of isocurvature perturbations in the two-field system of nonsingular Ekpyrotic 

scenarios. In particular isocurvature perturbations can be affected by the instability of the background near the 
bounce. In order to obtain the final spectra of adiabatic perturbations, we need to solve the coupled equations 
of adiabatic and entropy perturbations through the nonsingular bounce including radiation. It is important 
to investigate whether nearly scale-invariant spectra are obtained by conversion from isocurvature to adiabatic 
perturbations. 

Our analysis also applies to cyclic Universe models proposed in ref. [39,40] in which the bounce has been regular- 
ized by including higher-order corrections. Thus, our conclusions about the difficulty in obtaining a scale-invariant 
spectrum of fluctuations carry over to single field realizations of the cyclic scenario. In fact, we have done some 
simulations in the case of a simple negative potential V = m'^{(p — (pl) for |(^| < (j)c, and found that the solutions can 
be nonsingular so long as the higher-order effect dominates around the graceful exit. Note, however, that the spectra 
of density perturbations will be the same as in the Ekpyrotic scenario. 

Recently, a preprint has appeared [67] in which in the context of a brane world scenario a nonsingular bouncing 
cosmology is obtained by considering the motion of a D3 brane as a boundary of a five dimensional charged anti de 
Sitter black hole. In this model, computed in linear theory, the spectrum of gravitational wave fluctuations was shown 
not to be scale-invariant. This result supports the conclusions we have reached 

APPENDIX A: HEURISTIC DERIVATION OF THE SPECTRUM OF FLUCTUATIONS 

In this Appendix we give a heuristic derivation of the spectral index of cosmological perturbations in the PBB and 
Ekpyrotic scenarios. This analysis is based on two key assumptions. The first is the assumption that the amplitude 
of the fluctuations when they exit the Hubble radius during the phase of contraction (in the Einstein frame) is 
given by the Hubble constant. This assumption is reasonable assuming that the fluctuations are quantum vacuum 
perturbations which freeze when their wavelength crosses the Hubble radius. 

The second assumption is that the "physical magnitude" of the fluctuations remains unchanged while the wavelength 
of the fluctuation is larger than the Hubble radius. This assumption is much less obvious, although at flrst sight this 
assumption may seem obvious based on causality, namely the fact that microphysics cannot influence physics on scales 



'We are grateful to the Referee for pointing out this reference. 
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larger than the Hubble radius. However, in inflationary cosmology and in models with a contracting period such as the 
PBB and Ekpyrotic scenarios, the forward light cone (causal horizon) is much larger than the Hubble radius, and the 
spatial coherence of background fields over scales of the forward light cone can lead to nontrivial effects on fiuctuation 
modes on these scales, one of the most dramatic manifestations of this effect being the parametric amplification of 
super-Hubble (but sub-horizon) cosmological fluctuations during reheating in certain two-field inflationary models 
[68]- [74]. Furthermore, the term "physical magnitude" of cosmological fluctuations is not well determined. On super- 
Hubble scales, the magnitude of the density fluctuations depends sensitively on the coordinate system chosen. It is 
possible to choose coordinate-invariant (gauge-invariant) variables to describe the fluctuations, but there are many 
choices, and even in single field inflationary models many of these gauge-invariant fluctuation variables increase on 
super-Hubble scales (however, the increase between initial Hubble radius crossing during inflation at ti{k) and final 
Hubble radius crossing during the late time FRW cosmology at tf{k) is by a factor which only depends on the ratio 
of the equations of state at the two Hubble radius crossings). This increase is a self-gravitational effect. 

In spite of the above caveats, let us proceed with the heuristic discussion of the amplitude of density fiuctuations, 
applying it first to inflationary cosmology (exponential expansion to be specific). The quantity we wish to calculate 
is the mean square mass fluctuation on a scale k when the corresponding wavelength enters the Hubble radius at final 
Hubble radius crossing tf{k). This quantity, denoted \^ik,tf{k))\'^ is given by the power spectrum of fiuctuations 
[see eq. (3.20)], and its fc-dependence on the spectral index n is given by 



5M 



{k,tf{k)) 
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Based on the first assumption. 



SM 



{k,U{k)) 



H'^{U{k)) ~ const. . 



and using the second assumption we infer that 
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(6.4) 



(6.5) 



(6.6) 



and that hence the power spectrum is scale- invariant with an index n ~ 1. 

PBB cosmology is characterized (in the Einstein frame) by a scale factor which scales as 

a{t) ~ t^'^ , 

and thus 

H{t) = ^ ■ 

The condition of the initial Hubble radius crossing during the period of contraction 

ka-\ti{k)) = H{tiik)) , 

leads to 

ti{k) - , H{Uik)) ~ fc3/2 ^ 

and thus, applying our two basic assumptions as in the case of inflationary cosmology, to 

SM , 
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(6.7) 



(6.8) 



(6.9) 



(6.10) 



(6.11) 



which corresponds to a blue spectrum with index n = 4. 
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The analysis for Ekpyrotic cosmology is analogous. The only difference is that the value of p is different, < p <^ 1 
and hence 

fi(fc) ~ fc-Vd-p) , H{U{k))r^k^/^^-P\ (6.12) 

and thus, taking p = at the end, 

2 

~ , (6.13) 

which corresponds to a blue spectrum with index n = 3. 

Obviously, given the caveats discussed at the beginning of this Appendix, the results for PBB and Ekpyrotic 
cosmology cannot be trusted without a fully relativistic analysis. The growth rates of cosmological fluctuations are 
very different in expanding and contracting cosmologies, and thus even given that the above heuristic analysis works 
in the case of inflationary cosmology, this does not mean it has to work for PBB and Ekpyrotic cosmologies. However, 
the results of our paper are in agreement with those derived from the heuristic analysis. 
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APPENDIX B: ANALYTIC ESTIMATES FOR THE GRAVITATIONAL POTENTIAL 

Let us analyze the gravitational $ in more details. In the Einstein frame, the gravitational potential in the 
longitudinal gauge is expressed in terms of TZ^ in the absence of higher-order corrections [70,75]: 

*f = ^^f. (6.14) 

The gravitational potential in the string frame is related with the one in the Einstein frame as [76] 

^t = ^k+%^ = ^k-^Hk. (6.15) 

Making use of eqs. (2.10) and (4.10), we find that the curvature perturbation in the Einstein frame is exactly the same 
as that in the string frame (i.e, TZ^ = TZf). Therefore the gravitational potential in the Einstein frame is expressed 
in terms of Ttf : 

fl2 (ll^ -o/2 + oII^/2-M/4) 
4>f = — ^ ■■ -Tii . (6.16) 

Note that dots in eq. (6.16) denote the time-derivative with respect to ts- This is the equation which we solve 
numerically. 

Taking note of the relation, Hl{x) = H^-i{x) — {v/x)H^{x), one finds 



Aasz 



(6.17) 



where each sign corresponds to the case with 775 > and -qs < 0, respectively. When the evolution of z is given by 
z (X \ijsp , we have 

l-2r.T2hs|- = l-2^-27=l-27-|l-27| = |° tlZt]^;l ' (6-18) 



2(1 - 27) for 7 > 1/2 

This term completely vanishes during the contracting phase in the Ekyprotic cosmology with p < 1/3, since 7 is less 
than 1/2. In this case the gravitational potential in the Einstein frame can be expressed as 

This relation is used to estimate the spectrum of in Sec. IV. 
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